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We apply the technique of standard super symmetric factorization for any 
q factor ordering in the Wheeler-DeWitt (WDW) equation for barotropic 
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cosmological term. The resulting wave functions of the universe are exhibited 
as one-parameter families, which for particular values of the 7 parameter of 
the barotropic model are exactly found. 
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I. INTRODUCTION 



The study of eigenvalue problems associated with second-order differential operators 
t^. ; found a renewed impulse with the application of the factorization technique and its gener- 
alizations [1-4], which was fundamented on the inverse scattering method of Gelfand and 
^ ■ Levitan [5]. 

The use of the strictly isospectral scheme based on the general Riccati solution [1-4], 
has been applied from classical and quantum physics [1] to relativistic models [6]. This 
technique has been known since about a decade in one- dimensional supersymmetric quantum 
mechanics (SUSY-QM) and usually requires nodeless, normalizable states of Schrodinger 
equation. However, Pappademos, Sukhatme, and Pagnamenta [8] showed thah the strictly 
isospectral construction can also be performed on non-normalizable states. 

SUSY-QM may be considered an equivalent formulation of the Darboux transformation 
method, which is well-known in mathematics from the original paper of Darboux [7], book 
^ by Ince [9], and book by Matveev and Salle [10], where the method is widely used in the 

context of the soliton theory. An essential ingredient of the method is a particular choice of 
a transformation operator [5] in the form of a differential operator [1 1] which intertwines two 
hamiltonians and relates their eigenfunctions. When this approach is applied in quantum 
theory it allows one to generate a huge family of exactly solvable local potential starting 
with a given exactly solvable local potential [4] . 

In nonrelativistic one-dimensional supersymmetric quantum mechanics, the factorization 
technique was applied to the q = factor ordered WDW equation corresponding to the 
FRW cosmological models without matter field [12], where a one-parameter class of strictly 
isospectral cosmological FRW solutions was exactly found, representing the wave functions 
of the universe for that case. Also, in Ref. [13] a one-parameter family of closed, radiation- 
filled FRW quantum universe for any q factor order was found. 

In this work, we extend the application of the factorization technique to the case q 7^ 
factor ordering in the WDW equation for the barotropic FRW minisuperspace model, 
including a cosmological term. 
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The work is organized as follows. In section II we describe the supersymmetric fac- 
torization using the hamiltonian of the barotropic FRW cosmological model, including any 
factor ordering in the hamiltonian operator. In section III, we present the general solution 
to the Riccati equation obtained when supersymmeterizing this hamiltonian, as a strictly 
isospectral one-parameter family of both FRW cosmological potential and wave function of 
the model. Section IV is devoted to explicitely obtain the exact cosmological solutions of 
the WDW equation for some typical values of the 7 and k parameters in the literature. We 
give final remarks in section V. 

II. SUPERSYMMETRIC FACTORIZATION SCHEME 

We start with the hamiltonian that appears in the study of the barotropic FRW cosmo- 
logical model, with cosmological constant (Wheeler-DeWitt equation), 
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+ UAkA 2 + 48AA 4 - 384vrGM^A- 37+1 



|*) = 0, (1) 



where A is the scale factor, k is the curvature index of the universe ( k — 0, +1, —1 plane, 
close and open, respectively), A the cosmological term and 7 the parameter that describes 
the state equation of the fluid. 

In principle, the order ambiguity in equation (1) should be taken into account. This is 
quite a difficult problem to be treated in all its generality, since the hamiltonian operator in 
(1) must be written in a very general form in order to take into account all possible orders. 
For the semi-general Hartle-Hawking factor ordering [14], we have 1 

A-^^A-^A-»=A-f^- q A-»V (2) 



dA 2 dA dA V dA2 

where the real parameter q measures the ambiguity in the factor ordering. In this approach, 
the Wheeler-DeWitt equation can be written as follows 

d 2 ^ d# 

W = " A dA^ + q dA + V(A) * = °' (3) 

where 

V(A) = 144/tA 3 + 48AA 5 - 384ttGM 7 A" 37+2 . (4) 

Our interest here is to study the quantum solutions of equation (3) when using the 
supersymmetric factorization scheme. As we shall show, the solutions can be exactly found 
for particular values of the 7 parameter and the factor ordering parameter q. 

Consider the equation (modified WDW) 



1 There are other possibilities depending on different considerations on the operators, see refs. 
[15,16] 
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= -A- 2q ^| + q A-^ 2q ^ + A- 1 - 2q V(A)^ = 0. (5) 
dA 2 dA 

It is easy to show that the first order differential operators 

-4 + = -A- q ^- + W(A), (6) 

A~ = A-*±+W(A), (7) 

where W plays the role of a superpotencial function, factorize the hamiltonian (5) as 

H+ = A + A- . (8) 

The potential term V(A) is related to the superpotential function W(A) via the Ricatti 
equation 2 

V+(A, 7 ) = A^W 2 - A 1+q ^. (9) 



Making the transformation 



W 7 = -A- q ^, (10) 



u 7 



where the / means (9) is transformed into the original hamiltonian applied to the function 
u 7 , which means that the superpotential function is know once we have a solution to the 
original WDW equation. If we choose the ordinary factor ordering q = 0, we recover the 
relation found in [12]. 

In the supersymmetric factorization scheme, V_, the partner superpotential of V+, is 
obtained by performing the product 

FT = A~A + , H-f(A) = 0, (11) 

where f is the wave function related at the hamiltonian H~. Then, the isospectral potential 
to V+(A,7) is 

V_(A, 7) = A 1+2q W 2 + A 1+q w; = V+(A, 7) + 2A 1+q W 7 . (12) 

In addition, the solution f(A) is found using the fact that A + A~u^ = 0. Multipliying both 
sides of Eq. (12) by A~, we have that A~A + (A~u) = H (^4~u) = 0; then f(A) = A~ u. In 
this way, the structure for f(A) becomes 

f(A) = W 7 u 7 + A~ q u;. (13) 

Hence, knowing the superpotential function we can find the wave functions of the partner 
hamiltonian H~. However, this is not the most general solution, as will be shown in the 
following section. 



2 Hereafter we shall use W 7 to explicitely denote the superpotential functions depend on the actual 
choice of the 7 parameter, and we shall call V + to V in (5). 
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III. GENERAL SOLUTION 



The general solution to the Ricatti equation (12) is found when we propose the following 
scheme [1,4] 

V_ (A, 7) = A 1+2q W 2 + A 1+q W' = A 1+2q W 2 + A 1+q W 7 , (14) 
which by choosing 

W = W 7 + — , (15) 

leads to a Bernoulli equation for y 7 , 

y ;-2W 7 A q y 7 = A q , (16) 

whose solution is 

y 7 (A) = u~ 2 [I 7 + A] , (17) 

where I 7 (A) = x q u 2 dx. 

In this way, (15) can be written as 

W(A)=W 7 + I -^L (18) 

and the entire family of bosonic potentials can be built as 

V + (A, 7, A) = A 1+2q W 2 (A, 7, A) - A 1+q W'(A, 7, A), (19) 

V+(A, 7 ,A) = V_ -2A~ q W' (20) 

A 1+q u 7 u' A 1+2q < , x 

= V + (A, 7) - 4^-^ + 2^. (21) 

Finally, 



is the isospectral solution of the Schrodinger equation (3) for the new family potential (21), 
with the condition on the function g(A) = \J A(A + 1), though in the limit 

A — > ±00 g(A) — > A and u 7 — > u 7 . (23) 

This A parameter is included not for factorization reasons, because the wave functions 
in quantum cosmology are still nonnormalizable, but as decoherence parameter embodying 
a sort of quantum cosmological dissipation (or damping) distance. 

The WDW (3) equation has particular solutions for the 7 parameter and different uni- 
verses, which we shall explore in the following section. 
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IV. SOLUTION OF THE WDW EQUATION FOR PARTICULAR 7 



Quantum solutions can be readily found for particular choices of the 7 parameter. Here, 
we list some of them. 

1. 7 = — 1 corresponds to inflationary scenary. For this stadium (3) can be written 

Au" x - qu'_! + 144A 3 (m 2 A 2 - «)) u_i = 0, (24) 

where m 2 = -f + §7rGM_i. 

Now, we have three possible cases in our analysis: 
a) m 2 > 0, 

The differential equation for this subcase is 

Au" J - qu'_i + 144reA 3 (m 2 A 2 - k) u_i = 0, (25) 



which, after the consecutive substitutions v = m 2 A 2 — n, u_i = v 1 / 2 y(v), and z 
^jv 3//2 , yields an ordinary Bessel equation for q = 1, with solution 



u_i = (m 2 A 2 - k) 2 a Ji(z) +b J_i(z) 



(26) 



where z = [m 2 A 2 — k] , and a and b are superposition constants. 




2.b 



\ / 



FIG. 1. The isospectral wave function u_i (22) for m 2 = 4, ao = bo = 1. Curves corresponding 
to increasing A (solid line, A = 1, 11,61, 161,411) path u_i (dashed), eq. (26). 

b) m 2 < 0, 

The differential equation for this subcase is 

-Au" x + qu'.! + UAkA 3 (|m 2 |A 2 + k) u_i = 0, (27) 

being the solution the modified Bessel equation for q = 1 
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u_i = (|m 2 |A 2 + «;) 5 [aiIi(z)+biKi(z)] , (28) 

where z = t^jt [|?"n 2 |A 2 + k] 3 ^ 2 , and ai and bi are superposition constants, 
c) For m 2 = 0, and the differential equation for this situation is 

-Au'ii + qu'_i + 144/cA 3 u_i = 0, (29) 

which, for k — 1 has as solutions the modified Bessel functions of order v = i±a 

u_! = A 2 - [A L,(6A 2 ) + B K,(6A 2 )] , (30) 

where A and B are superposition constants, whereas for k = — 1, the solutions 
become the ordinary Bessel functions 

u_i = A 2u [AxJ^GA 2 ) + BiYj,(6A 2 )J , (31) 

where Ai and Bi are superposition constants. 

2. Dust era, 7 = 0, and plane universe, k = 0: 

By use of the transformations z = |>/3A A 3 and uo = e~ z / 2 w(z), we find the hyperge- 
ometric differential equation for w(z) 

d 2 w . dw „ 

z-TV + a-z— -nw = 0, (32) 
az z dz 

where n = ^ — 1 ^ r ™° and a = Therefore, the independent solutions are [17] 

wi = iFi(n,a;z), (33) 
w 2 = z 1 - Q 1 F 1 (n-a + l,2-a;z), (34) 

where 1F1 is the degenerate hypergeometric function. In this way, the solution for this 
case become 

u = e- z / 2 [A 2 w 1 (z) + B 2 w 2 (z)], (35) 

where A 2 and B 2 are superposition constants. 

3. Stiff fluid, 7 = 1 and plane universe 

The WDW equation for this stadium is written as 

Au'l - qu; + 48 (-AA 5 + SttGMiA -1 ) m = 0, (36) 
whose solutions become 
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) 2 - 384ttGM 



(37) 



where Z M is a generic Bessel function with real or imaginary order /i [18]. 

In this case, the exact expression of the solutions will depend on the signs of the 
cosmological constant and the parameter fi. 

i) /x real and A > 0, the function become the modified Bessel functions, either I M 
or K M , depending on the boundary conditions. 

ii) ii real and A < 0, the functions Z M turn into the ordinary Bessel function, either J M 
or Y M ; depending on the boundary conditions. 

iii) A > and \i pure imaginary, the functions Z M become the modified Bessel functions 
of pure imaginary order [18], either I M or K M , depending on the boundary conditions. 

iv) A < and fj, pure imaginary, the functions Z M turn into the ordinary Bessel 
functions of pure imaginary order [18], either J M or Y M , depending on the boundary 
conditions. 



SUSY-QM has allowed us to show that new classes of exact solutions are found for the 
WDW equation, extending the class of exactly solvable spectral problem for the Schrodinger 
like hamiltonian in one space dimension (in quantum cosmology the WDW equation plays 
this role). In this work, we found isospectral cosmological potential and a one-parameter 
family of wavefunctions of the universe for barotropic FRW models with a cosmological 
constant, including the factor ordering problem. We were able to find exact particular 
solutions for different choices of the 7 parameter in the WDW equation (3) and different 
universes. The parameter A looks like a decoherence parameter embodying a sort of quantum 
cosmological dissipation (or damping) distance (see Fig. (1)). 
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